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Abstract 

In the annual reports other aspects of materials processing of high temperature materials, 
such as Heat Transfer in Fiber Formation, have been discussed. In this Final Report Axial 
Segregation in Unsteady Diffusion Dominated Directional Solidification of a Binary Alloy 
with a Planar Solid-Liquid Interface Profile is considered. The analysis is an extension of 
the work of Smith ct al. [l] to account for an arbitrary initial axial concentration profile. 
This situation is important in the experimental studies of growth of crystals in a micro gravity 
environment in which the growth is from melts formed by re-melting crystals that have a 
known axial impurity profile. Three situations are discussed. These are the initial transient, 
a re-solidified crystal, and multipass solidification. 
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1 Introduction 

In this article axial segregation in unsteady diffusion dominated solidification of a binary 
alloy is considered when an orfciirurj/ initial concentration profile is given. It is an extension 
of the classical study of Smith et al. [1] who give the solution for dopant concentration 
distributions during initial and final transients in crystals grown from melts in which impurity 
distribution is initially uniform- The analysis should be of general interest to crystal growers, 
but the particular needs of experimenters who carry out crystal growing experiments in a 
microgravity environment were kept in mind. In their experiments crystals are grown from 
a melt that has been formed by re-melting a crystal that was grown by some method on 
earth. Hence the initial axial impurity distribution in the melt is not uniform, but a specified 
function of the axial coordinate. 

In coder to carry out the analysis by analytical methods, a number of assumtions were 
made. The major ones include neglecting convection and changes in the volume as the 
material melts. The first is closely approximated in space, and the second is an assumed 
property of the material. In addition rcmclting between successive passes was assumed to 
be instantanous and thus diffusion in the melt was neglected during the remelting process. 
How well this can be achieved in practice depends on the furnace design and the smallness 
of the diffusion coefficient. The final transient was also ignored, since it contributes only 
to the local concentration in the far end of the ampoule. The influence of phase boundary 
curvature was neglected in this study. Its influence under steady growth has been most 
recently reported by Korpela et al. [2]. 

From the above list, it is dear that in order to carry the analytical study a number of 
aasuptions are needed, assumptions which could have been relaxed, in a purely computational 
study. The advantage of carrying the analysis analytically is the possibility of obtaining 
explicit solutions which are easily programmed and from which special results for various 
materials can be gleaned. It is with this aim that this study was undertaken. 

2 Formulation 

Consider the directional solidification of a melt in an infinite cylindrical ampoule. The 
periphery is insulated and impermeable, and the temperatures of the ends are such that 
solidification takes place. The origin of coordinates is fixed to the mean position of the 
solidification front. If the coordinates move with the solidification rate, in these coordinates 
the solidification front is stationary and the molt moves toward it at a constant speed V. 

For semiconductors the coefficient of mass diffusion D in the solid is several orders of 
magnitude smaller than it is in the fluid phase. Thus solute diffusion In the solid can be 
ignored. By introducing the non-dimensional variables 
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the concentration field of the solute is governed by 

9C ac &C 

dt dx ~ dx* w 

and initial, boundary, and interface conditions are 

C(x,0) */(*), C(oo,t)-l, |£(0,t) + (l-*)C(0,i)=0 (3) 

where the distribution coefficient k = C,(0)/C(0, t) has been introduced. It is the ratio of 
the slopes of the liquidus and solidus lines in an ideal phase diagram in which these lines are 
taken to be straight. The initial condition /(x) is a given function that tends to unity far 
from the solidification front. 


3 Solution 


The mathematical problem formed by Eqns. (2) and (3) for unsteady diffusion equation is 
solved by employing a Green’s function, which satisfies 

dG_dG PG 

dt dx~ dx a W 

C(x,0)=*(x-0, — (0,t) + (l-k)G?(0,t)-0, 0(00,0=0 (5) 

The solution obtained by Laplace transform technique. The Grcen’B function is 

C(i| (,t) = (i - t) - (I - 1) vsl 


Although the operator for the Green’s function is not self-adjoint, the solution for the con- 
centration field can be expressed, as 6hown in Stalcgold [3], in term of the original Green’s 
function by proper identification of the variables x and (. Details of the solution method are 
given in Kim [4]. 

In terms of the Green’s function the concentration field can be written as 

e(*,g-jfG(*ie,‘)/<0* m 

This is a convenient form not only for obtaining analytical solutions for simple initial concen- 
tration fields, but also for numerical evaluation for quits general initial distributions /(x). 
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4 Results 


4.1 Initial Transient 


When the initial melt is uniform in composition no that /(x) = 1. integration of the Green’s 
function gives 


CM- l G(x|(,t)iif = * er fc[_ij(i_()]_i e ,fc[_i^( I + t)| 

This result is not new, for it was first obtained by Smith ct al. [l]. It is given here as a 
special ease of Eqn. (7). In dimensional variables the concentration distribution is 

£i dr 1 ~ 1 + - v<, » - wm 

- (1 - »)Vf)] (») 

The concentration at x* = Vf multiplied by h is frozen into the crystal because D,j D ^ 
1. Thua the impurity concentration in the crystal in dimensional terms, obtained by putting 
X* — Q and then setting f * x'/V, is 


c.M-kCi 0,t') = ^ 


[1 


+ erf(^^)] + (2it - l)e-< 1 -*l fcV *'^erfc[(21t - 1) 



(io) 

This result for C, is shown in Figure 1 for various values of k, and it Is the same as that first 
given by Smith at al. [1]. It is repeated here for the sales of completeness. 


4.2 Re-melting and Solidification 

In space experiments crystals are grown by re- solidifying a melt that has been formed from 
re- melted crystals. Hence for them the initial condition /(x) =* C t [x) is that obtained from 
the initial transient solution (10), which in dimensionless terma is 

/(*) -;[[! + ««V?)] + (» - l)e- (1 -‘ )b, *rfc[(» - 1)^]] (11) 

The concentration in the re-solidified crystal is given by Eqn. (7) with thin as the initial 
condition. Thus the impurity concentration in the crystal be expressed as 


C,(t) - *C(U,t) = ^f° 


(i _ *) e -o-t>*4i* 


eric| i7? - - i)v7 l + 
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(1 + Qrf(y|)] + (2k - erfc[(2fc - 1)^|]| (12) 

The axial distribution in the solidified crystal is now obtained by simply replacing t by x in 
this expression and it is given by 

-If [<5 - Q'-"- 1 *"* -(!-*)/=] + 

[1 + ert(<J\)} + (2k - erfcpfc - 1)^|1 if 


(13) 

The solution is in the form of a quadrature and is plotted in. Figure 2 for different valu& of 
k. 

If k ■* 1/2, the solution simplifies greatly for three of the terms making up the integrand 
vanish and the remaining terms can be integrated in closed form. By splitting the integral 
in 


m 


di 


(14) 


into two parts, the first one is 

jf = yfnx ef ^1 + erf(^)] 

Tb evaluate the second, the error function is represented by its series 

„ , 2 ^ (-l)"a**+i 

V*£* «!(2n+l) 

and a term by term integration is carried out. This yields 

f ei_fi Wf* ~ * ) f e4 ~ £ ^ 

The integral on the right can be written as 

f «*-* +|:|; f) 

and U(a\ b; c) is the confluent hypergemuetric function of the second kind [5). 
The concentration distribution in the crystal can thus be written as 


(IS) 


(16) 


(17) 


( 18 ) 


*■<*>- i 


!+ «rf (\j-) 
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or in dimensional terms as 


CW} 1 

C. 4 


l+Kf (vSr> 


I d* f (-D W )* 1 "** 1 . 3 n 3 1 tV 

rr £o 2 an+3 (2n + l)n! 11 + 2^2 + 4’ 2’ 4 D* ' 3 °' 


This solution is plotted in Figure 3 and it together with Figures 1 and 2 shows the effect 
of a second pass on the reduction of the impurity concentration. 


4.3 Multipass Solidification 

Often in industrial practice there is a need to purify crystals by sone refinement. Although 
only a region of finite length is ranched, the theory in this section serves as a qualitative 
guide to the purification achieved by multiple meltings and subsequent solidification*. 

After the first pass concentration in the solid can be written as 

C.to * JT G(0\& x)/(Odf • (21) 

Since the initial condition /(£) is the previous solid concentration, for the n-th pass the 
concentration can be expressed as 

c!”M - jTc( oif.xjej- 1 )©*. (22) 

Evaluating the integral numerically gives 

c^M = £) c(ofc,* i )c<- ,) KfH, (23) 

J=1 

where wj's are the weigths of a given quadrature rule. In matrix notation this can be written 

(24) 

where the matrix elements are Ay = Gy • wj (no summation convention on indices is implied 
here). Using this result recursively gives 

CS’ - K,0 S’ (25) 

whan CS> is the initial concentration distribution. This is a handy formula, which allows the 
concentration for the n-th pass of solidification to be obtained directly by matr ix multipli- 
cations. Wide availability of computational tods, such as Matlab, for matrix calculations 
makes generation of results particularly simple. The results shown in Figure 4 are for up to 
50 re-eolldifieatiozv. 
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5 Conclusions 

By finding the Green’s function for unsteady concentration equation, analytical results were 
obtained for the concentration distribution in a crystal, which initially has a uniform impurity 
distribution, and for the distribution in a re-solidified crystal, with an initial concentration 
distribution set by the first pass. This advances the theory over that worked out by Smith 
et al. [1] by allowing arbitrary initial conditions to be considered. The results will serve 
the needs of experimentalists who carry out experiments in a micro- gravity environment. 
Since the solution (7) is in terms of a Green’s function, the n-th pass of re-solidification can 
be evaluated by numerical integration and a few matrix multiplications. This mak& the 
analysis of multiple passes particularly simple to program and thus useful to practitioners. 
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List of Symbols 


Symbol 

Representation 


concentration of impurity 

Co 

concentration of impurity in the far field 

D 

diffusion coefficient 

G 

Green's function 

k 

partition coefficient 

t 

time 

U 

confluent hypergeometric function 

V 

velocity 

X 

spatial coordinate 

r 

Gamma function 


location of delta function 


List of Figures 

Pig. 1 . Ax ia l impurity on initial transient solidification for various values of k. 

Fig. 2. Axial impurity on re-solidification for various values of k. 

Fig. 3. Effect of re-solidification on impurity concentration of crystal for k =* 0.5. 

Fig. 4. Effect of number of solidification on impurity concentration of crystal for k = 0.5. 
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Figuie 1. The influence of the initial transient on the axial impurity concentration for various 
values of k. 
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Figure 2: Axial impurity concentration after re-sojidification for various values of k 
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Figure 3: Effect of re-solidification on impurity concentration for a crystal with k — 0 . 5 . 
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